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Using lattice simulations of quenched QCD we estimate the 
finite size effects present when a gluon plasma equilibrates in 
a slab geometry, i.e., finite width but large transverse dimen- 
sions. Significant differences are observed in the free energy 
density for the slab when compared with bulk behavior. A 
small shift in the critical temperature is also seen. The free 
energy required to liberate heavy quarks relative to bulk is 
measured using Polyakov loops; the additional free energy re- 
quired is on the order of 30 - 40 MeV at 2 - 3 T c . 



I. INTRODUCTION 

The formation of a quark-gluon plasma in a central 
heavy ion collision is generally assumed to take place 
in a coin-shaped region roughly 1 fermi in width, with 
radius comparable to the radii of the colliding nuclei, 
which is to say several fermi. While lattice gauge theory 
has given us information about bulk thermodynamic be- 
havior, finite size effects have up to now been studied 
using simplified, phenomenological models. We study 
via lattice gauge theory simulations the behavior of a 
gluon plasma restricted to a slab geometry, with the lon- 
gitudinal width much smaller than the transverse direc- 
tions. This inner region is heated to temperatures above 
the bulk deconfinement temperature, surrounded by an 
outer region which is kept at a temperature below the 
deconfinement temperature, providing confining bound- 
ary conditions for the inner region. Details of this work 
are given in |jj . 

Measurements of the equilibrium surface tension cuo of 
pure SU(3) lattice gauge theory (quenched lattice QCD) 
show that the dimensionless ratio cto /Tj 3 is small. For the 
case of an TV 3 x 4 lattice, a /T c 3 w 0.0292(22). § A sim- 
ple estimate of surface tension effects on the transition 
temperature can be obtained from a simplified model in 
which only volume and surface terms appear, as in the 
bag model. Q Based on these simple considerations, fi- 
nite size effects due to the surface tension should be small. 

Other contributions to finite size effects come from 
a variety of sources. In the case of systems with non- 
abelian symmetries, global color invariance produces an 
additional finite volume effect which will not be consid- 
ered here. j| H M In general, finite size effects lead 
to a rounding of the transition. This can be taken 
into account in the bag model by a Maxwell construc- 



tion, leading to mixed phases and a broadened critical 
region. A recent treatment for the quark gluon plasma 
can be found in || . We have attempted to avoid these fi- 
nite volume effects by making the transverse dimensions 
large. 



II. METHODOLOGY 

In lattice calculations, finite temperature is introduced 
by the choice of N t , the extent of the lattice in the (Eu- 
clidean) temporal direction. The relation of physical 
temperature T to N t and the lattice spacing a is sim- 
ply T — l/N t a. The lattice spacing a implicitly depends 
on the gauge coupling (3 = 6/g 2 in a way determined by 
the renormalization group equations. To lowest order in 
perturbation theory, the relation is given by 
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where A^ is renormalization group invariant and the 
renormalization group coefficients 6q an d b\ are given by 
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In analyzing our data, we used the renormalization 
group results given in reference ||, which are deter- 
mined directly from lattice simulations, and contain 
non-perturbative information about the renormalization 
group flow. 

By allowing the coupling constant [3 to vary with spa- 
tial location, a spatially dependent temperature can be 
introduced. We have chosen the temperature interface to 
be sharp, in such a way that the lattice is divided into two 
spatial regions, one hotter and one colder. The quenched 
approximation simplifies the role of the cold region, be- 
cause below T c , the dominant excitation at low energies 
is the scalar glueball. Temperatures near T c are smaller 
than glueball masses by about a factor of four, so glue- 
balls play no essential role in the thermodynamics, and 
the pressure in the hadronic phase is essentially zero. We 
thus expect the slab thermodynamics to be largely insen- 
sitive to the precise temperature of the region outside the 
slab, as long as it is sufficiently low. In full QCD, this in- 
sensitivity to the outer temperature would not hold, due 
to pions. Note that the role of boundary conditions here 
is quite different from that in bubble nucleation. In that 
case, both (3 in and [3 out are taken to be near T c . [fill fTO] 



1 



The slab is given a fixed lattice width, w = 6a, rather 
than of fixed physical width. Since N t — 4, wT is fixed 
at 3/2. At higher temperatures, the width in physical 
units is somewhat smaller than the longitudinal size of 
the plasma formation region expected in heavy ion colli- 
sions. While the use of equilibrium statistical mechanics 
to study gluon plasma properties during the early stages 
of plasma formation may appear suspect, a simple esti- 
mate using the Bjorken model jl^] shows that when a 
coin-shaped region of width 1 fermi has expanded to 1.5 
fermi, the variation in temperature is only from 0.8T at 
the center of the coin to T at its edges. 

The free energy density / for the slab is obtained using 
the standard method |13[ of integrating the lattice action 
with respect to (3. We use a convenient convention for 
the sign of / that is opposite the usual one. In the bulk 
case, / is then identical to the pressure p. 
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where S = (\/N)ReTr U p . As in the bulk case, it is 
necessary to subtract the zero-temperature expectation 
value from the finite temperature expectation value, in 
this case using the same pair of values. In general, 
quantum field theories with boundaries develop diver- 
gences that are not present in infinite volume or with pe- 
riodic boundary conditions. Such divergences would re- 
quire additional boundary counterterms. Symanzik 
has shown to all orders in perturbation theory that in the 
case of 4> 4 with so-called Schrodinger functional bound- 
ary conditions that the theory is finite in perturbation 
theory after adding all possible boundary counterterms 
of dimension d < 3 consistent with the symmetries of the 
theory. It is generally believed that this result applies 
as well to all renormalizable field theories and general 
boundary conditions, but a proof is lacking. Luscher et 
al. [[[5) have shown for gauge theories that at one loop 
no new divergences are introduced by Schrodinger func- 
tional boundary conditions. This is consistent with the 
non-existence of gauge-invariant local fields of dimension 
< 3 in pure Yang-Mills theory. 

In order to take advantage of the data on bulk thermo- 
dynamics provided by the Bielefeld group || , we worked 
consistently with lattices of overall size 16 3 x 4. The 
values used for each subtraction come from 16 4 lattices 
with identical values of /3j„ and (3 ou t ■ The value of f3 ou t 
was held fixed at 5.6 while varied from 5.6 to 6.3. 
For comparison, the bulk transition for iVf = 4 occurs at 
(3 c (N t = 4,N S = 16) = 5.6908(2) f3 c {N t = 4, N s = oo) = 
5.6925(5). § 



III. FREE ENERGY OF GLUONS 

Figure 1 shows the free energy density f/T 4 versus 
T/T c compared with the bulk pressure. The free en- 
ergy in the slab is lower than the bulk value by almost 



a factor of two at 2T C . It appears that the slab value is 
slowly approaching the bulk value, but other behaviors 
are also possible. Calculations of the finite-temperature 
contribution to the Casimir effect for a free Bose field 
contained between two plates show that f/T 4 has a non- 
trivial dependence on the dimensionless combination wT. 
Jl6| [jl7| It is natural to ask if the corrections to the free 
energy seen here can be accounted for by the conven- 
tional Casimir effect. A straightforward calculation of 
the free energy of a non-interacting gluon gas confined to 
a slab shows an increase in the free energy density over 
the bulk value by a factor of about 1.63 at wT = 3/2. 
The Casimir effect cannot explain the reduction of the 
free energy observed in our simulations. 
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FIG. 1. Slab free energy compared with bulk pressure. 

A consistency check was performed on the surface ef- 
fects. JXs| The free energy density was calculated for a 
system at /3j„ = 6.0 by performing simulations with (3 in 
fixed at 6.0 and (3 out varying from 5.6 to 6.0. Combin- 
ing these results with the bulk data of reference || cre- 
ates a path equivalent to varying /3, n while holding (3 ou t 
fixed. For (3 in = 6.0 and (3 out = 5.6, this gives f/T 4 = 
0.65 ± 0.04, to be compared with f/T 4 = 0.69 ± 0.03 for 
the direct calculation. The major source of systematic 
error lies with the choice of boundary conditions for the 
slab, here set by Pout- We have estimated the effects of 
varying j3 ou t by performing simulations at = 6.2 and 
Pout — 5.5 on 16 3 x 4 and 16 4 lattices. These results sug- 
gest that lowering f3 ou t from 5.6 to 5.5 reduces the free 
energy by roughly 10 percent at (3i n =6.2. 



IV. SURFACE TENSION 

We define an effective surface tension a(w, T) by 

f = p-2a(w,T)/w (4) 

where the notation a(w,T) recognizes that the surface 
tension a does depend on the width of the slab and the 
internal and external temperatures. The factor of 2 oc- 
curs because the slab has two faces. In the limits where 
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T approaches T c and w goes to infinity, this quantity ap- 
proaches aQ. Figure 2 shows a(w, T)/T 3 versus T/T c for 
wT = 3/2; representative error bars are shown. 
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divergence can be eliminated when comparing bulk ex- 
pectation values to those in finite geometries. We define 



AF Q (x) = -T In [P s i ab (x)/P bu , 



(7) 



as the excess free energy required to liberate a heavy 
quark in the slab geometry relative to bulk quark matter 
at the same temperature. This technique can also be 
used in, e.j.^a spherical geometry, which is relevant for 
nucleation. || |(J 

In figure 3, we show the expectation value for the 
Polyakov loop versus z measured in lattice units for sev- 
eral values of f3. Each curve is normalized by dividing 
the values of the Polyakov loop by the bulk expectation 
value at the corresponding value of /3. Error bars are 
shown only for even values of z. 



FIG. 2. Non-equilibrium surface tension versus T. 

The value at /3 = 5.7, 0.056 ± 0.002 is higher than the 
value a = 0.0292 ± 0.0022 given in reference |§. We at- 
tribute this to two effects: in our case (3 ou t is fixed at 5.6, 
whereas for equilibrium measurements it is extrapolated 
to P c , and our finite value of the width w also acts to 
increase a(w, T) over the equilibrium value as measured 
in simulations at large w. Away from the bulk critical 
point, a/T 3 rises quickly to a peak at about 1.4T C , and 
then falls slowly as T increases. A large non-equilibrium 
surface tension has also been observed in measurements 
of the equilibrium surface tension, where these effects 
were obstacles to obtaining ctQ. 



V. FREE ENERGY OF QUARKS 



The Polyakov loop defined by 
P{x) = (l/N c )Tr Texp 



l/T 

i I Ao [x, t) dr 
o 



(5) 



is the order parameter for the deconfinement transition 
in pure (quenched) gauge theories. In the case of SU(N), 
there is a global Z(N) symmetry which ensures at low 
temperature that the expectation value (TrP) is 0. At 
sufficiently high temperatures, this symmetry is sponta- 
neously broken. The expectation value of the Polyakov 
loop can also be interpreted in terms of the free energy 
of an isolated, infinitely heavy quark Fq: 



(P(x))=exp[-F Q (x)/T] 



(6) 



In the low-temperature confined phase, Fq is taken to 
be infinite, whereas in the high temperature phase it is 
finite. Direct extraction of Fq from computer simula- 
tions is problematic, because the expectation value has a 
multiplicative, /3-dependent ultraviolet divergence. This 




10 12 14 16 



FIG. 3. Polyakov loop versus z. 

It is clear that a significant change occurs between 
13 = 5.8 (T = 1.23 T c ) and [3 = 5.85 (T = 1.36 T c ). 
For larger values of (3, AFq diminishes to a value of ap- 
proximately 30 — 40 MeV in the middle of the slab. 
In table 1, we list (3, T, slab width in fermis, width of 
the core in fermis, and AFq in MeV for representative 
values. The width of the core is calculated by interpolat- 
ing the Polyakov loop profiles and determining the region 
where the slab expectation value is greater than 80% of 
the bulk value. All conversions to physical units are per- 
formed by taking the string tension a to be (425 MeV) 2 , 
which implies T c = 254 MeV. § 

TABLE I. /3, T, slab width in fermis, width of the core in 
fermis, and AFq in MeV for representative values. 
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T (MeV) 


w (fm.) 


W core (fm.) 


AFq (MeV) 


5.8 


313 


0.95 





169 


5.85 


346 


0.86 


0.40 


57 


6.0 


455 


0.65 


0.48 


41 


6.2 


624 


0.47 


0.40 


31 
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VI. CONCLUSIONS 

There are significant deviations in the slab geome- 
try from bulk behavior and the ideal gas law, arising 
from a strong non-equilibrium surface-tension. This non- 
equilibrium surface tension can be an order of magnitude 
greater than the equilibrium value. Surface tension ef- 
fects also produce a mild elevation of the apparent crit- 
ical temperature. Measurement of Polyakov loop expec- 
tation values relative to bulk shows that the suppression 
of heavy quark production due to the slab geometry is 
small. 

There are good reasons to call our result an estimate 
rather than a calculation. Although lattice gauge theory 
simulations of bulk behavior can be made arbitrarily ac- 
curate in principle, in this case there is some uncertainty 
in the precise exterior boundary conditions appropriate, 
and indeed in the applicability of equilibrium thermody- 
namics at this early stage of quark-gluon plasma forma- 
tion. However, this seems like the best estimate available 
now, and further refinements are possible. 

We have not yet explored the nature of the phase tran- 
sition, which will require some care. One interesting pos- 
sibility is that the order of the transition might change as 
the width changes. The deconfinemcnt transition in bulk 
quenched finite temperature QCD is in the universality 
class of the three-dimensional three-state Potts model, 
which has a first-order phase transition. As the width 
of the slab becomes commensurate with the correlation 
length near T c , the phase transition should cross over to 
the universality class of the two-dimensional three-state 
Potts model. The two-dimensional three-state Potts 
model has a second-order phase transition, so it is possi- 
ble that the order of the transition may change. |2Q] The 
correlation length at the bulk transition is known to be 
large pH , so it is likely that the transverse correlation 
length in the gluonic sector is much larger in the slab 
geometry than in bulk, even if crossover does not take 
place. 

The system studied here has some interesting features 
amenable to theoretical analysis. In the outer region, 
the Polyakov loop will decay away from the interface as 
exp[— ar/T out ] , where a is the string tension. In the inner 
region, the Debye screening length sets the scale for the 
Polyakov loop. We are currently working on a theoretical 
model of this system based on perturbation theory which 
has these features. 
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